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?H ■ ABSTRACT 

c/3 i Dynamical models for 17 early-type galaxies in the Coma cluster are presented. The 

C3 . galaxy sample consists of flattened, rotating as well as non-rotating early-types includ- 

ing cD and SO galaxies with luminosities between Mb = —18.79 and Mb = —22.56. 
I \ Kinematical long-slit observations cover at least the major and minor axis and extend 

■ to 1 — 4 r c g ■ Axisymmetric Schwarzschild models are used to derive stellar mass-to- 
y—{ ' light ratios and dark halo parameters. In every galaxy the best fit with dark matter 
Q\ , matches the data better than the best fit without. The statistical significance is over 

■ 95 percent for 8 galaxies, around 90 percent for 5 galaxies and for four galaxies it is 
' not significant. For the highly significant cases systematic deviations between observed 
, and modelled kinematics are clearly seen; for the remaining galaxies differences are 

■ more statistical in nature. Best-fit models contain 10-50 percent dark matter inside 
I the half-light radius. The central dark matter density is at least one order of magni- 
. tude lower than the luminous mass density, independent of the assumed dark matter 

density profile. The central phase-space density of dark matter is often orders of mag- 
• i-H , nitude lower than in the luminous component, especially when the halo core radius is 

■ large. The orbital system of the stars along the major-axis is slightly dominated by 
' radial motions. Some galaxies show tangential anisotropy along the minor-axis, which 

is correlated with the minor-axis Gauss-Hermite coefficient H4. Changing the balance 
between data-fit and regularisation constraints does not change the reconstructed mass 
structure significantly: model anisotropies tend to strengthen if the weight on regular- 
isation is reduced, but the general property of a galaxy to be radially or tangentially 
anisotropic, respectively, does not change. This paper is aimed to set the basis for 
a subsequent detailed analysis of luminous and dark matter scaling relations, orbital 
dynamics and stellar populations. 
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1 INTRODUCTION 



Elliptical galaxies are numerous among the brightest galax- 
ies and they harbour a signific ant fraction of the present- 
day s t ellar mass in t he universe (|Fukugita. Hogan &i Peebles! 
1998; R enzinil [20061 ). Key parameters for the understand- 
ing of elliptical galaxy formation and evolution are, among 
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others, the central dark matter density, the scaling ra- 
dius of dark matter, the stellar mass-to- light ratio and 
the distribution of stellar orbits. While the concentration 
of the dark matter halo puts constrai n ts on the assembly 
epoch llNavarro. Frenk fc White! Il996l ; Ijing fc SutJ |2000| ; 
IWechsler et al.l 2002f ). the orbital state co ntains imprints of 
the assemb ly me chanism of ellipticals fe.g.lvan Albadalll982l; 
HernauistJ fl99^. Il993|; IWeil fc HernauistJ Il996l; iDubinskil 



1998l ; lNaab fc Burkerti2003l ; |jesseit. Naab fc Burkertll2 005) 
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Information about elliptical galaxy masses are in 
principle offered through various channels. The analy- 
sis of X-ray halo temperatures, the kinematics of occa- 
sional gas discs and galaxy-galaxy lensing provide ev- 
idence for extended dark matter halos around early- 
type galaxies (e.g. iBertola et al] fl993l: IPizzella et al] 



19971; lLoewenstein'fc White! 1 19991; lOosterloo et al] |2002|: 



Hoekstra et aljfeoO^Fukazawa et al]|2006l;lHumphrey et al] 
20061 ; iKleinheinrich et al] l200rj ; iMandelbaum et al.l bo06). 



These methods do not constrain the inner halo-profiles 
strongly, however. At non-local redshifts strong lensing con- 
figurations allow a detailed reconstruction of the mass en- 
close d inside, say, r c s (e.g. iKeetonl l200ll ; iKoopmans et al] 
2006). None of the above mentioned observational channels 
is sensitive to dynamical galaxy parameters, such as the dis- 
tribution of stellar orbits. 

Dynamical modelling of stellar kinematics has the 
unique advantage that it allows reconstruction of both the 
mass structure and the orbital state of a galaxy. High- 
quality observations of the line-of-sight velocity distribu- 
tions (LOSVDs) out to several r c g are needed for this pur- 
pose. To overcome the problems of measuring absorption 
line kinematics in the faint outskirts of ellipticals, discrete 
kinematical tracers such as planetary nebulae or globular 
clusters can be used to additiona l ly constrain the mass dis- 



tribution (e.g. Saglia et al] |2000| ; iRomanowskv et al] 120031 ; 



iPierce et alj|2006h 

Since stars in galaxies behave collisionlessly to first or- 
der, the distribution of stellar orbits is not known a priori 
and very general dynamical methods are required to probe 
all the degrees of freedom in the orbital system. So far only 
one large sample of 21 round, non-rotating giant ellipticals 
has been probed for dark m atter considering a t least the full 
range of spherical models (jKronawitter et al]|2000f ). These 
models predict circular velocity curves constant to about 10 
per cent and equal luminous and dark matter somewhere in- 
side 1 — 3 r e ff . Reconstructed halos of these models are ~ 25 
times denser than in comparably bright spir als, which indi- 
cates a ~ 3 times higher formation redshift l|Gerhard et al] 
l200ll ). Not all apparently round objects need to be intrinsi- 
cally spherical; some may be face-on flattened systems. 

Apparently flattened ellipticals have not yet been ad- 
dressed in much generality. Primarily, because axisymmet- 
ric modelling is required to account for intrinsic flattening, 
inclination effects and rotation. Fully general axisymmetric 
models involve three integrals of motion, one of which - the 
non-classical so-called third integral - is not given explicitly 
in most astrophysically relevant potentials. Only recently, 
sophisticated numerical met hods such as Schwarz schild's or- 
bit superposition technique fSchwarzschild 1979) have pro- 
vided fully general models involving all relevant integrals of 
motion. Dynamical studies of samples of elliptical galaxies 
using this tech nique are, however, based on kinematical data 
inside r < r eB i|Gebhardt et alj|2003l ; ICapperlari et al]|2006l ) 
and dark matter is not considered. 

The present paper is part of a project aimed to analyse 
the luminous and dark matter distributions as well as the 
orbital structure in a sample of flattened Coma ellipticals. 
The data for this project has been collected over the last 
years and consists of ground-based as well as (archival and 
new) HST imaging and measurements of line-of-sight ve- 
locity distributions (LOSVDs) along various position angles 



out to 1 - 4r».fT ([Mehlert et al.l [2OO0I ; IWegner et ail |2002| ; 
ICorsini et al. 1 12007^ 7 The implementation of our modelling 
machinery, which is an adv anced version of the axisym - 
met ric Schwarzschild code of iRichstone fc Tremaine] (Il988l ) 



and iGebhardt et al l (|2000|) has been described in detail in 
iThomas et al. I i|2004 I 20051 ) . In the present paper we sur- 
vey the models of the whole sample. This sets the basis 
for subsequent investigations of luminous and dark matter 
scaling relations and stellar populations in elliptical galaxies 
(Thomas et al. 2007a, in preparation). 

In Sec.[5]the observations are summarised and the mod- 
elling is outlined in Sees. [3] and U The mass structure of our 
models and the orbital anisotropies are described in Sees. [5] 
and [()] respectively. Phase-space distribution functions for 
luminous and dark matter are the subject of Sees. [7] and [5] 
We discuss the influence of regularisation on our results in 
Sec. [5] The paper closes with a short discussion and sum- 
mary in Sec. 1101 A detailed comparison of models and data 
for each galaxy can be found in App. IA1 



2 SUMMARY OF OBSERVATIONS 

The Coma sample consists of seventeen early-type galax- 
ies: two cD galaxies, nine ordinary giant ellipticals and six 
lenticulars or galaxies of intermediate type. They cover the 
luminosity interval —20.30 < Mb < —22.56, typical for lu- 
minous giant ellipticals/cDs. One single fainter galaxy with 
Mb = —18.8 is also included (cf. Tab. [1] magnitudes are 
from Hyperleda for a distance of d = 100 Mpc to Coma; 
this corresponds to Ho = 69km/s/Mpc). Effective radii are 
mostly between 3'.'3 < r c s < 18'.'4. Only the four bright- 
est galaxies have fo rmally very large r e g ~ 30" — 70" (cf. 



Tab.[T] r c ff are from I Mehlert et al] l|200ol ) and based on de- 
Vaucouleurs fits). All galaxies share the same distance and 
the spatial resolution in the photometric as well as the kine- 
matical observations is roughly comparable for all galaxies. 

The photometric input for the modelling is constructed 
as a composite of ground-based (outer parts) and HST imag- 
ing (inner parts). The two surface brightness profiles /j gr d 
and /ihst are joined by shifting the HST profile according 
to the average (/i gr d — ^hst) over a region where both data 
sets overlap and seeing effects are negligible (~ 5" — 16"). 
The shift {fJ, gv d — Mhst) is usually well defined (cf. Tab.[T|. 

The kinematic data to be fit by the models come from 
long-slit observations along at least two position angles: the 
apparent major and minor axis, respectively. Kinematical 
parameters from different sides of a galaxy are averaged. 
As error-bars we use the maximum from the two sides or 
half of the scatter between them, whatever is larger. This 
assumes that uncertainties in the observations are mostly 
systematic (see also the discussion in Sec. I4.1|l . In case of 
pure statistical errors and an exactly axisymmetric system 
this would overestimate the errors by a factor y/2. Thus, we 

are conservative. 

T he data are de s cribed in full detail in [Mehlert et al] 

|2000l ). IWegner et al] i|2002l ) and lCorsini et al] l|2007l ). Basic 
parameters of the photometric and kinematic data set are 
summarised in Tab. [T] 

Three galaxies deserve further comments: 
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Tabl e 1. Summary of ob servational data. (1, 2) galaxy id (GMP from [Godwin. Metcalfe &i Peach 1983]); (3) morphological type 
(from lMehlert et aljlioool) ; (4, 5) HST and ground-based photometry (L97 = HST/WFPC2 R-band data, Principal Investigator: 
John Lucey, Proposal ID: 5997; H98 = HST/WFPC2 R-band data, Principal Investigator: William Harris, Proposal ID: 6104; W07 
= HST/ WFPC2 R-band data , Principal Investigator: Gary Weg ner, Proposal ID: 1 0884; MOO = Kron-Cousins il^-band photom- 
etry from lMehlert et aljlipool: J94 = Gunn r photometry from Ijcirgensen & Franx 199 41); (6) absolute B- band magnitude (from 
Hyperleda; Ho = 69km/s/Mpc); (7, 8) effective radius r c g and ellipticity e e at r e g from lMehlert et al.l d2000] ) ; (9) rms(/i gr( j — Mhst) 
between shifted HST surface brightness /^hst and corresponding ground-based /x gr( j; (10-13) radius of the outermost kinematic 
data point along various slit positions: maj/min/dia = position angle of 0°/90°/45° relative to major axis; off = parallel to major 
axis (in case of GMP5568: two offset-slits). The offsets are quoted in the captions of Figs. IAT1 - IA17I 



GMP5568/NGC4816: GMP5568 has been observed 
along four position angles. In addition to major and minor- 
axis spectra, two observations were made with the slits par- 
allel to the major-axis: one with an offset of r e ff /4, the other 
with r cff /20. 

GMP0144/NGC4957: The velocity dispersion peak of 
GMP0144 is significantly off the photometric centre. Fur- 
thermore, GMP0144 is the only galaxy in our sample that 
exhibits a significant isophotal twist towards the centre. 
Thus, GMP0144 is likely triaxial near its centre. To reduce 
the influence of potentially non-axisymmetric regions on our 
modelling, kinematic measurements inside r < 4" are omit- 
ted. 

GMP5975/NGC4807: Dynamical models for GMP5975, 
based on major a nd minor- axis kinem atics, have already 
been presented in iThomas et al.l (I2005T ). There, a striking 
depopulation of retrograde orbits was found. To check its 
significance we also determined kinematics along a diagonal 
slit. Here we present new models that include this additional 
kinematical data. Both, the mass structure and the distri- 
bution of stellar orbits did not change significantly. 



3 DYNAMICAL MODELLING 

We model the kinematic and photometric obs ervations with 
Schw arzschild's orbit superposition method (|Schwarzschildl 
19791). Details about our implementation are given in 
Thomas et al.1 l|2004 120051 ). Basic steps of the method are 
briefly recalled in this section. 



3.1 Deprojection and inclination 

The surface photometry is deprojected to obtain the 3d lu- 
m inosity distributi on v for each galaxy (using the program 
of iMagorrianl [199 9) . We consider radial profiles of surface- 
bri ghtness, ellipticity and isopho tal shape parameters and 
06 l|Bender fc M ollcnhoff 119871 ) for the deprojection^ (cf. 
App.[A}. For each galaxy, we probe three different inclina- 
tions in the deprojection, and subsequent dynamical mod- 
elling, respectively: (1) i — 90° (edge-on), where the de- 
projection is intrinsically least flattened; (2) a minimum in- 
clination that is found by requiring the deprojection to be 
intrinsically as flattened as an E7 galaxy; (3) an interme- 
diate inclination, for which the deprojection looks like an 
E5 galaxy from the side. This inclination scheme emerges 
as a compromise between limited computation time on the 
one side and the strategy to get the most conservative es- 
timate of uncertainties on intrinsic properties on the other. 
In many galaxies the inclination is only poorly constrained 
(cf. Sec. 14. 2[) . In other words, when we quote 68 percent 
confidence uncertainties on intrinsic properties below this 
includes in many cases models from all three probed incli- 
nations, including those connected with the rather extreme 
intrinsic E5 and E7 shapes. 

In case of GMP0756, GMP1176, GMP1990 and 
GMP2417 only the edge-on orientation is considered. These 
galaxies are all highly flattened. In addition, they appear 
either discy (e.g. GMP1176) or have thin dust features 

1 In case of GMP1176 isophotal shape parameters up to a\i 
are used to repre sent the isophotes appropriately (see also 
ICorsini et"alll2007l) . 
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(GMP1990, GMP2417; cf. ICorsini et~ai1 12007ft , implying 

that they are seen close to edge-on. 

The dynamical modelling of GMP5975 in lThomas et al.l 

(|2005l ) revealed that only models at i = 90° were within 
the one sigma confidence region. It was also found that the 
deprojection of GMP5975 becomes implausibly boxy in the 
outer parts, if it is assumed that the galaxy is significantly 
inclined. We therefore reanalysed the extended data set for 
GMP5975 only with i = 90°. 



3.2 Mass model 

With the luminosity density v given, a trial mass density 
distribution can be defined by 



p = T v + pdm- 



(1) 



The stellar mass-to-light ratio T is assumed constant 
throughout the galaxy. Concerning the dark matter density 
Pdm we probe the two following parametric prescriptions. 
First ly the NFW-distribution jNavarro. Frenk fc White! 
1 19961 ) 



PNFw(r, r s ,c) oc 



(r/r a )(l + r/r a )^ 



(2) 



where r a = r2oo/c is a scaling radius, r2oo is a mea- 
sure of the virial radius and c is the concentration of 
the halo. Simulations predict the two halo parameters r s 
and c to be correlated, such that the distribution ((2| can 
be read as a one-parameter f amily of dark matter halos 
l|Navarro. Frenk fc Wliitelll996h . To be explicit, we use 



r*. oc io( A - logc )/ s 



(3) 



with A = 1-05 a nd B = 0.15 (|Navarro. Frenk fc White1l996l ; 
iRix et al]|l997h . We consider spherical as well as flattened 
NFW halos, where the latter are derived from equation ((2} 
by the substitution r — » r y/ cos 2 (#) + sin 2 (#) /q 2 (q is the 
constant flattening of the isodensity contours, $ is the lati- 
tude in spherical coordinates). 

The second halo family used is the logarithmic potential 



PlogM oc v\ 



3rj + r 
(r 2 + r 2 ) 



(4) 



that gives rise to an asymptotically constant ci rcular velocit y 
v c and a flat central density core inside r < r c |Binnevll98 ll ) . 

In the gravitational potential generated by the mass 
distribution Jl] ) we comput e typic ally about 18000 orbits as 
described in lThomas et alj l|2004 ). 



and / being the phase-space distribution function (DF) of 
the model. In Schwarzschild models - by construction - the 
DF is constant along individual orbits. The corresponding 
phase-space density along orbit i is the ratio 



f — — 

H = V t 



(7) 



of the total amount of light Wi on the orbit (the so-called 
orbital weight) and the orbital phase-space volume Vi. The 
Wi that solve equation ((5J are obtained iteratively: starting 
with a low a = 10~ 10 we solve equation (J3J for a fixed set 
of oti, using the orbital weights obtained at as initial 

guess for the solution at a%. 

The x 2 -te-mi in equation © quantifies deviations be- 
tween model and data. We do not include the photometry 
in the x 2 > but treat the deprojected luminosity distribu- 
tion as a boundary condition for the solution of equation 
J5J. To fit the measured LOSVDs, which are parameterised 
in t erms of the G a uss-Hermite parameters v, a , H3 and 
H 4 i| Gerhard! 1 19931 ; Ivan der Marel fc Franxl Il993l ) we pro- 
ceed as follows: the Gauss-Hermite parameters are used to 
generate binned data LOSVDs £^ ar Errors are propagated 
via Monte-Carlo simulations. These data LOSVDs and the 
corresponding model quantities C^ od are used to get the 
Xlosvd of equation ©: 



Xlosvd 



EE 

3=1 fe=l 



r? - c- 

mod dat 



AC 



(8) 



The above sum includes all Nc data points and each LOSVD 
is represented by N ve \ bins in projected (line-of-sight) veloc- 
ity. 

With the orbital weights u?j determined, the dynamical 
state of the model is completely specified, i.e. the phase- 
space distribution function is known (cf. equation [7]). In the 
course of this paper we will not only consider the DF, but 
also the orbital anisotropy. It can be quantified by the so- 
called anisotropy parameters 



a» = 1 - -i 

at 

(meridional anisotropy) and 



0* 



1 



(9) 



(10) 



(azimuthal anisotropy) . Internal velocity dispersions a in the 
above equations are computed in spherical coordinates r, -Q 
and ifi, oriented such that <p is the azimuth in the equatorial 
plane and 1? is the latitude. 



3.3 Orbit superposition 

The final orbit superposition model is constructed 
according to the maxim um entropy technique of 
iRichstone fc Tremaind (| 198Sf ) . It consists in solving 

S = S -a xlosvd -> max, (5) 
with S denoting the Boltzmann entropy 

S = - f f]n{f) d 3 rd 3 ^=-^^ln(^) (6) 



3.4 Regularisation 

The parameter a in equation J5} controls the relative weight 
of data-fit and entropy maximisation. The higher a the bet- 
ter the fit and the larger the noise in the derived distribution 
function, or orbital weights, respectively. Ideally, regularisa- 
tion has to be optimised case-by-case for each galaxy. This 
holds in principle for both, the value of the regularisation 
parameter a, as well as for the functional form of S. Firstly, 
because the spatial resolution and/or coverage as well as 
the signal-to-noise of the observations vary from galaxy to 
galaxy and regularisation should be adapted to that. This 
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NFW 


-1.37 


19.1 


90 



Table 2. Summary of modelling results. (1) galaxy id (cf. Tab. [TJ; (2-3) best-fit stellar Tsc [Mq/Lq] (ij^-band) and achieved 
goodness-of-fit Xsc (P er data point) without dark matter; (4-7) as (2-3), but for logarithmic halos with parameters r c [kpc] and v c 
[km/s]; (8-11) same as (2-3), but for NFW halos with concentration c and flattening q; (12-13) best-fit halo-profile with significance 
A xLlo = (Xnfw ~~ X LO g) x ^data; (14) evidence for dark matter Ax^, M = (Xg C ~ XmiJ x ^data; (15) inclination of best-fit with 
minimum and maximum in the 68 percent confidence region of calculated models (where no range is quoted, only edge-on models 
were calculated). 



primarily concerns the choice of a. Secondly, different galax- 
ies have different intrinsic structures. Specifically, the degree 
to which the entropy of a stellar system is maximised may 
vary in phase-space. Consider, for example, a cold disc in- 
side a hot spheroid. The disc has low entropy and to fit its 
rotation, a needs to be large (models with the maximum 
entropy according to equation [6] have no rotation). On the 
other hand, the spheroid-dominated region in phase-space 
can have higher entropy and using a large a in the fit ampli- 
fies the noise in the corresponding parts of the phase-space 
distribution function (DF). The dilemma as to the choice of 
q in such cases could be solved by adjusting the function S 
appropriately. 

For the Coma galaxy modelling we use the same regu- 
larisation for all galaxies: a — 0.02 and the entropy of equa- 
tion ([6]). The value for a has been obtained by Monte-Carlo 
simulations of isot ropic rotator test ga laxies with realistic, 
noisy mock data (|Thomas et al.ll2005l ). Applying it to the 
whole sample is motivated by the similar spatial coverage 
and resolution of all our Coma observations (cf. Sec. [2}. Fur- 
thermore, since it has been obtained from fitting isotropic 
rotators it has proven to be sufficiently large to fit non- 
maximum entropy, rotating galaxies. It might be slightly 
too large for non-rotating galaxies. Thus we expect models 
of non-rotating galaxies to possibly be noisier than those 
of rotating galaxies, but we do not expect that our imposed 
regularisation is too restrictive. In any case, we will explicitly 
investigate the dependency of model results on the choice of 
q in Sec. [9] 



3.5 Best-fit model and uncertainties 

To obtain the best-fit mass model for a given a we calculate 
orbit models as described in Sees. 13.21 and 13.31 for various 
combinations of the relevant parameters: (r c ,v c ,T) in case 



of LOG-halos, (c, T, q) in case of NFW-halos or just T for 
models without dark matter, respectively. Logarithmic halos 
are probed on a grid with Ar c ~ r c ff/2 and Av c = 50 km/s 
(in some cases the grid is refined around the location of the 
lowest \ 2 ). Typically we explore N r x N v ~ 100 halos. The 
analogous numbers for the NFW halos read Ac = 2.5, N c ~ 
12 and N q = 2 (q G {0.7, 1.0}). The step-size AT for the 
mass-to- light ratio equals 10-20 percent of the best-fit T c i yn , 
independent of the halo type. Around the best-fit model the 
resolution in T is doubled, resulting in Nx ~ 6 — 10 models 
with different mass-to-light ratios for each halo. This sums 
up to about 600 — 1000 models with logarithmic halos and 
120 - 240 models with NFW halos. The final number of 
models is up to a factor of three larger, depending on the 
number Ni ^ 3 of probed inclinations i. The total number 
of models per galaxy is around 1000 — 3000. 

Among these models we determine the best-fit accord- 
ing to the lowest Xgh 5 defined as 



N C 



Xg 



3=1 
JJ 3,mod 



A "dat 



JJ 3,dat 



+ 



A °dat 



+ 



+ 



H 3 



H 3 



(11) 



Here, v 3 d&t is the rotation according to the Gauss-Hermite 
parameterisation of the LOSVDs (other parameters analo- 
gously). A detailed discussion ab out the differen c e betw een 
Xgh an d Xlosvd can be found in lThomas et al.l (|2005l ). 

Confidence intervals on model quantities are set by the 
corresponding minimum and maximum values obtained over 
all models within A^gh = 1.1 from the minimum Xgh- The 
value of Axgh = IT i s slightly more conservative than the 
classical A^gh = 1-0 an d has been derived by means of 
Monte-Carlo simulations (|Thomas et al.ll2005l ). 
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Figure 1. Comparison of confidence intervals with (shaded) and without (dotted) rescaling the \ 2 (details in the text). Individual 
panels show circular velocity D c ; rc , mass densities p (stellar: blue/dotted, dark matter: black/dashed and total: red/solid), dark 
matter fraction Mdm/M, Xqh (dotted without rescaling, solid with rescaling) and velocity anisotropies j3$ and fi v along major 
and minor-axis, respectively. 



As a byproduct of the iterative technique to solve equa- 
tion lf5)l. we get - for each set of (r c , v c , T, i), (c, q, T, i) and 
(T,i), respectively - orbit models for about N a m 50 dif- 
ferent values of the regularisation parameter. This allows us 
to derive a best-fit model for each 014 and to explore the 
dependency of best-fit model parameters on a (cf. Sec. O. 



4 MODELLING RESULTS 

Modelling results are summarised in Tab. [2] In the remain- 
der of this Section we collect some general notes on these 
results. 

4.1 Goodness-of-fit 

The goodness-of-fit obtained under the different assump- 
tions about the overall mass distribution are given in 
columns (3), (7) and (11) of Tab. H Thereby 

xlc = min{xGH(T, i)/iV data }, (12) 

Xlog = mimtxcH^c, v c , T, i)/7V data } (13) 

and 

Xnfw = min{xGH(c, q, T, i)/iVdat a } (14) 

are minimised over all relevant mass parameters. Differences 
between models with and without dark matter are further 
discussed in Sec. 15.11 Here we only refer to the fact that our 
models are able to reproduce the observations with a Xgh 
per data point which is in many cases significantly smaller 
than unity. The largest deviations between model and data 
occur for the SO GMP0756, possibly related to the low H4 
along the offset-axis, which are not followed by our models. 
Fits to some galaxies are as good as Xmin i$) 0.1, where 



Xmin = min{xLOG, Xnfw.Xsc} (15) 

describes the overall minimum of Xgh f° r a given galaxy. In 
many cases where Xmin is particularly low, error bars of the 
observations are much larger than the point-to-point scat- 
ter of the data points. This concerns GMP5279, GMP2921, 
GMP4928, GMP5568 and GMP3958, where the observa- 
tional errors are likely overestimated (see also the fits in 
App.|A|. In some other systems, like for example GMP0144, 
the error bars used in the modelling are rather large, because 
they also include side-to-side variations of the kinematics, 
which are often also larger than the point-to-point scatter 
on a given side of the galaxy. Thus, both effects might partly 
explain the low Xmin °f these galaxies. 

Very low Xmin raise the question whether confidence in- 
tervals of model properties calculated as described in Sec. 13. 51 
(and shown in Figs.[4][5][6]and[7]below) are overestimated. In 
cases where the observational errors are obviously too large 
it is reasonable to rescale the m until y^ iT1 ~ 1 . In fa ct, this 
has been done for GMP5975 in lThomas et all (|2005l ). where 
error bars were scaled such that Xmin ~ 0.7. This value was 
determined from Monte-Carlo simulations of isotropic rota- 
tor models. To quantify the effect of rescaling, Fig. [1] ex- 
emplifies confidence intervals for one galaxy of the sample 
(GMP4928) once without rescaling the Xgh an d once with 
rescaling all observational error bars to Xmin — 0.7. As it 
can be seen, uncertainty regions shrink a lot after rescaling. 

Globally rescaling the error-bars is not appropriate 
for all galaxies, however. Errors along the major-axis of 
GMP3792, for example, might be overestimated, but those 
along the minor-axis seem not. Moreover, in a case like 
GMP1750 slight minor-axis rotation - which cannot be fit 
with axisymmetric models - adds a constant to Xmin- Just 
rescaling to Xmin ~ 0.7 in all galaxies would introduce an ar- 
tificial dependency of uncertainty regions on minor-axis ro- 
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Figure 2. Confidence levels A%q h versus T (normalised to the best-fit T^yn). Solid: logarithmic halos; dashed: NFW halos; 
dotted: mass follows light; horizontal dashed: 90 percent confidence limit. Where a dotted line is missing the self-consistent 
case is ruled out with more than 95 percent confidence (exception: GMP3329, where the best- fit Tgrj ~ 1.7T^ yn is outside the 
plotted region). From top- left to bottom-right galaxies are arranged in order of decreasing total mass inside r e g . 



tation, which is not appropriate. In order to treat all galaxies 
of the sample homogeneously we do not rescale the Xgh> but 
give the most conservative error-bars for our models. The 
corresponding confidence intervals may be interpreted as the 
maximal uncertainty on derived model quantities, while the 
shaded regions of Fig. Q] may be interpreted as lower limits 
for these uncertainties. 

Fig. [2] shows the dependency of 



Axgh(Y) = xgh(T) - min (xgh), 



(16) 



where minxcm = Xmin x ^data (cf. equation^ and Xgh( t ) 
is minimised over all NFW- fits, logarithmic-halo fits and 
self-consistent fits with the given T. For all but one galaxy, 
we find a clear minimum in Xgh CO- The exceptional case, 
GMP3329, is peculiar in many respects: (1) It is among the 
brightest galaxies of the sample and has a very large r e g. 
The data only cover the region inside r < r e ff/2. (2) The 
surface-brightness profile shows a break near 0.2 — 0.3 r e ff 
(cf. upper panel of Fig. IA1|) . (3) At about the same pro- 
jected distance from the centre the velocity dispersion dips 
and rises again at larger radii (cf. lower panel of Fig. IA1|) . 
The poor constraints on the mass-to-light ratio in this sys- 
tem could be related to the poor data coverage. It might 
also be that GMP3329 is actually composed of two subcom- 
ponents. If these have different mass-to-light ratios T, then 
the XGH _curve mav have two corresponding local minima 
and the poor data coverage may smooth out these into a 



flat plateau. Finally, our modelling of GMP3329 may suf- 
fer from the Coma core being possibly not in dynamical 
equilibrium, as ind icated by the kinem atics of intra-cluster 
planetary nebulae (|Gerhard et al.ll2007T ). 



4.2 Model inclinations 

Most of the best-fit models are edge-on according to the last 
column of Tab. [5] This is surprising at first sight because if 
galaxies are oriented randomly then one would expect only 
about 6 out of 17 objects to have inclinations larger than 
i > 70°. Omitting the five systems for which only edge-on 
models were calculated (GMP0756, GMP1176, GMP1990, 
GMP2417 and GMP5975; cf. Sec. E3J and taking into ac- 
count the uncertainties quoted in Tab.[2]there are 3 galaxies 
out of 11 where inclinations i < 70° are ruled out by our 
modelling (at the 68 percent confidence level). This is in 
good agreement with the expectation for random inclina- 
tions. Nevertheless, we now discuss a little more whether 
our modelling might be subject to a slight inclination bias. 

First, one possible issue is that using the same regu- 
larisation for all galaxies might introduce a subtle bias to- 
wards edge-on configurations. Consider a rotating system: 
the lower the assumed inclination of the model the larger its 
intrinsic rotation needs to be in order to match the observa- 
tions after projection. Thus, the system will be dynamically 
colder and its entropy will be lower. Turning the argument 
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Figure 3. Top: histogram of apparent short-to-long axis ratios 
at r e fj. Bottom: intrinsic best-fit short-to-long axis ratio (b/a) 
(averaged over r/r e ff £ [0.5,2.5]). Black/solid: whole sample; 
red/dashed: without SOs. 



around: usage of a constant a enforces the same weight on 
entropy in the inclined model as in the edge-on model. Since 
the inclined model has to have lower entropy, however, its fit 
may be less good. This might drive models of rotating galax- 
ies towards i — 90°. We do not expect this to affect conclu- 
sions drawn from our modelling results strongly, because, as 
it has been argued in Sec. 13.11 error bars on intrinsic prop- 
erties include in many cases models at different inclinations, 
even extreme cases. But it might drive the best-fit model to 
occur preferentially around i — 90°. 

Second, for face-on galaxies noise in the kinematics 
may be a source of bias towards edge-on models as well 
|Thomas et al.ll2007bl ). It can cause rotation measurements 
v 7^ even for exactly face-on, axisymmetric galaxies. An 
edge-on model can in principle fit these » / 0, whereas face- 
on models necessarily obey v = 0. Thus, everything else 
fitting equally well, the contribution of noise in v to the \ 2 
would be smaller in edge-on than in face-on models. Since 
we have no clear candidate face-on galaxy in our sample we 
do not expect this issue to be relevant to the Coma sample, 
however. 

The third thing to note is that galaxies for our sam- 
ple may be not at random inclinations. The sample is 
designed to complement earlier studies on round galaxies 
and we explicitly selected flattened, rotating ellipticals and 
SOs to be fitted. The distribution of apparent axis-ratios 
(b/a) e = 1 — e e (with e e from Tab.[T]) is shown in the upper 
panel of Fig. [3] The sample exhibits a tail of highly flat- 
tened systems, which lacks, for example, in the distribution 
of bright galaxies with de-Vaucouleurs pro files in the Sloan 
Digital Sky Survey l| Vincent fe Rvdenl2005h . This tail is pro- 
duced by the SO galaxies in our sample and clearly shows 
that the sample as a whole is biased towards flattened sys- 
tems. The ellipticity distribution of those galaxies that are 



classified as ellipticals in Tab. [T] (dashed line in Fig. [3} is 
still shifted slightly towar ds higher ellipticit i es wit h respect 
to the bright ellipticals of lVincent fc Rvdenl (|2005h . In com- 
bination with the lack of round objects in our sample, this 
indicates that even our 11 ordinary ellipticals are slightly 
biased, but the sample is too small for a definite appraisal. 

Fourth, even assuming our sample galaxies are at ran- 
dom inclinations and that the first two just discussed points 
are irrelevant (regularisation and noise) and that inclina- 
tions can be reconstructed uniquely from ideal data with 
ideal models then we still could be faced with a slight bias 
in our models: as it has been described in Sec. l3.ll we do not 
probe a fine grid in inclinations but look for extreme cases. 
Our models provide for each galaxy only the choice between 
edge-on or intrinsically E5/E7, respectively. Since an intrin- 
sic E5/E7 shape is a rather extreme assumption this might 
drive the modelling towards the edge-on option as well. 

The distribution of intrinsic axis ratios of the Coma 
models is shown in the lower panel of Fig. [3] It peaks 
at b/a — 0.8, consistent with deprojections of the fre- 
quency function of ellipti c al galaxy apparent Bat tenings 
i|Tremblav fc Merrittl [l996l : I Vincent fc Rvdenl 120051 ). Com- 
pared with these studies, the distribution in the lower panel 
of Fig.0has relatively more galaxies on the flatter side of the 
peak and relatively fewer galaxies on the rounder side. Now, 
if the modelling would be subject to a strong bias towards 
i — > 90° then we would expect the opposite: an overestima- 
tion of intrinsically roundish galaxies. Thus, the lower panel 
of Fig. [3] argues against a strong modelling bias towards 
high inclinations. However, the argument is not conclusive, 
because the sample itself maybe biased against apparently 
round galaxies. This could partly compensate for an incli- 
nation bias in the modelling. 

In conclusion, although there might be a slight inclina- 
tion bias in the modelling and/or the sample galaxies, Fig. [3] 
reveals that either this bias is not very strong, or that mod- 
elling and sample biases roughly counterbalance each other. 



5 LUMINOUS AND DARK MATTER 
DISTRIBUTION 

Now we discuss the distribution of luminous and dark matter 
in the Coma models. 



5.1 Does mass follow light? 

According to Tab. [51 the best-fit model of each galaxy 
contains a dark matter halo. Column (14) of the ta- 
ble states that eight galaxies have at least a two sigma 
detection of a dark matter halo (GMP0282, GMP0756, 
GMP1176, GMP1750, GMP2440, GMP4928, GMP5568 and 
GMP5975). The best fitting models with and without dark 
matter, respectively, are compared to the kinematic data in 
App. [X] From this comparison it follows that models with- 
out a halo obviously fail to reproduce the kinematic data for 
the above mentioned galaxies. The evidence for dark mat- 
ter thereby comes mostly from the innermost and outermost 
kinematic data points: without dark matter, the energy of 
the models is too low, when compared to the data at large 
radii and too high, when compared to the central data (for 
example illustrated by the dispersion profile of GMP5975). 
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Figure 4. Compilation of circular velocity curves. Thick/red: luminous+dark (68 per cent confidence region shaded); dot- 
ted/blue: luminous only; dashed: dark matter only; vertical dotted lines: boundaries of kinematical data; arrows denote the 
effective radius r cf r. From top-left to bottom-right galaxies are arranged in order of decreasing total mass inside r ef j. 



The reason for the differences at small radii is likely that 
part of the missing outer mass in models without a halo is 
compensated for by a larger mass-to- light ratio (cf. Sec. l5.5p . 
This, in turn, causes an increase of the central mass and 
central velocity dispersion, respectively. In GMP1750 the 
dispersion profile without dark matter fits systematically 
worse than the one with dark matter at all radii. Concerning 
GMP5568, the dispersion along one of the offset-slits is par- 
ticularly large, larger than in all other slits. It is not entirely 
clear if these large dispersions are real. If not, then they er- 
roneously increase the evidence for dark matter. However, 
because the error-bars of the corresponding data points are 
very large, these dispersions are not the dominant driver for 
the dark halo detection in GMP5568. 

In the rest of the sample the detection of dark matter - 
if at all - is more of statistical nature. Models with and with- 
out dark matter for GMP0144 and GMP2921 differ in a sim- 
ilar fashion as those of GMP1750. For GMP3510, GMP3958 
and GMP5279 small differences between models with and 
without dark matter can be seen at the last kinematic data 
points, but the formal significance for dark matter is less 
than 90 percent. We believe that the statistical significance 
for dark matter in these five cases is underestimated due to 
our very conservative error estimates. 

In the four remaining objects GMP1990, GMP2417, 
GMP3329 and GMP3792 the evidence for dark matter is 
generally low. Poor evidence for dark matter in GMP3329 
maybe related to the overall poor constraints that the mea- 
sured kinematics put on its mass-to-light ratio (cf. Sec. l4.1[) . 



GMP1990 is consistent with the assumption that mass fol- 
lows light. 

Our sample thus roughly divides into three categories: 
(1) galaxies that are clearly inconsistent with a constant 
mass-to- light ratio (8 galaxies out of 17). (2) Cases where 
models with and without a dark halo differ systematically, 
but where the formal evidence for a dark halo is less than 
two sigma (5 galaxies). In these cases we expect that our 
very conservative error bars lead to an underestimation of 
the dark matter detection. (3) Systems in which the evidence 
for dark matter is generally weak (4 galaxies). 

Models and data of some galaxies with a clear dark 
halo detection still differ systematically in the outer parts 
(e.g. GMP0756, GMP1176 and GMP5975). Decreasing the 
weight on regularisation reduces these differences. However, 
according to the discussion of Sec. 13.41 we do not expect 
that we have significantly over-regularised our models. Even 
in case we would have, the derived halos of these systems do 
not depend much on the choice of the regularisation param- 
eter, such that conclusions upon the masses of these galaxies 
are robust (cf. Sec. [5}. It might be possible that differences 
between models and data are related to changes in the stellar 
population, that other our adopted halo profiles are not ap- 
propriate for these systems, or that the corresponding outer 
regions are not in equilibrium or not axisymmetric. We plan 
further investigations of these topics for future publications. 
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Figure 5. Top: spherically averaged mass densities. Red (solid): total mass; blue (dotted): luminous mass; black (dashed): 
dark matter with errors (shaded). Bottom: corresponding spherically averaged dark matter fractions with 68 percent confidence 
regions. Vertical dotted lines and arrows as in Fig. [4] 
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5.2 Circular velocity curves 

Fig. [4] shows the best-fit circular velocity curves for the 
Coma galaxies. The shapes of the curves vary from cases 
with two local extrema (e.g. GMP5279) to a case of mono- 
tonic increase (GMP3958). The sample provides examples 
of rising as well as falling outer circular-velocity curves. 
Flattened, rotating galaxies have fairly flat circular velocity 
curves beyond the central rise (less than 10 percent variation 
up to the last kinematic data point in GMP0282, GMP0756, 
GMP2417, GMP3510, GMP1176 and GMP5975). 



5.3 Mass densities and dark matter fractions 

Spherically averaged density profiles of all Coma galaxies 
are surveyed in the top panel of Fig. The luminosity 
distribution of most galaxies shows a power-law core that 
smoothly joins with the outer light-distribution. Towards 
the most luminous galaxies the central slope of the luminos- 
ity distribution flattens out (GMP3329 and GMP2921). The 
inner breaks in the light profiles of GMP1990 (r ~ 2kpc), 
GMP2417 (r « 1 kpc) and GMP2440 (r « 0.3 kpc) originate 
from prominent dust features. 

The central regions are dominated by luminous matter. 
Halo densities in the centre are at least one order of mag- 
nitude lower than stellar mass densities - independently of 
the halo profile being either of the logarithmic or of the 
NFW type. The radius where dark and luminous mass den- 
sities equalise is inside the kinematically sampled region of 
each galaxy. In some galaxies the transition from the lumi- 
nous inner parts to the dark matter dominated outskirts is 
very smooth (for example GMP0756, GMP1176, GMP5975). 
The corresponding dark halo components are relatively con- 
centrated (NFW halos) and the circular velocity curves are 
fairly flat. In other galaxies the transition is marked by a 
break in the total mass profile and a dip in v c i rc (for exam- 
ple GMP0282). We will come back to these different circular 
velocity curve shapes in Thomas et al. (2007a, in prepara- 
tion). 

Dark matter fractions of the best-fit Coma models are 
shown in the bottom panel of Fig. [5] In most galaxies 10 to 
50 percent of the mass inside r e g is dark. 



5.4 NFW or logarithmic halos? 

The evidence for or against logarithmic and NFW halos is 
summarised in column (13) of Tab. [2] The majority of best- 
fit models (13 out of 17) is obtained with logarithmic halos. 
However, the significance for one or the other halo profile 
providing the better fit is in most cases low. As already 
implied by the approximate flatness of the circular velocity 
curves, the overall effect of the halo component is to keep 
the outer logarithmic slope of the total mass density around 
—2 (i.e. the case of an exactly constant v c irc). This can be 
achieved either with logarithmic halos (asymptotically) or 
with suitably scaled NFW halos (over a finite radial range 
around the scaling radius). Differences in the profiles' inner 
slopes seem to play a minor role, perhaps because the inner 
mass profile turns out to closely follow the light profile. If 
elliptical galaxy circular velocity curves are roughly flat over 
a very extended radial range, then NFW fits will break down 



at some point. With the data at hand no clear decision in 
favour of one of the two profiles can be made. 

Most of the best-fit NFW-models are obtained with a 
flattened halo. Since we cannot significantly discriminate be- 
tween NFW and LOG-halos, firm statements about the flat- 
tening of the halos are not possible. 

5.5 The fraction of mass that follows the light 

From Tab.[5]it can be taken that mass-to-light ratios Tsc of 
self-consistent models are on average (17±10) percent larger 
than those of models with a dark matter halo. Concerning 
the difference between logarithmic and NFW halos, the best- 
fit Tnfw is generally equal or lower than the corresponding 
Tlog- 



6 VELOCITY ANISOTROPY 

Having explored the mass structure of the models we next 
focus on their dynamics. In this section we will consider 
the velocity anisotropics defined in equations (|9} and (|10[) . 
respectively. 

The discussion will be restricted to the minor-axis and 
major-axis bins of the Schwarzschild models, respectively. 
Each galaxy is covered by kinematical observations along at 
least these axes and the internal orbital structure is best 
constrained there. 

6.1 The polar region 

Fig.[6]surveys velocity anisotropy profiles along the intrinsic 
minor-axis of the Coma galaxy sample. According to axial 
symmetry fi$ = f3 v holds directly on the symmetry axis. 
Hence, the upper and lower panel of Fig. [6] are overall very 
similar. The minor-axis bins of the models, however, form 
a cone with opening angle Ai? = 25° around the z-axis. 
Thus, they include regions off the symmetry axis, where the 
equivalence between azimuthal and meridional dispersions 
does not hold, such that fig and fi v in Fig.[6]are not identical. 

In Fig. [6] the very central anisotropies should not be 
regarded as reliable. Firstly, because the central bins are 
affected from incomplete orbit samp ling, resulting in ar - 
tificialty large azimuthal dispersions l|Thomas et alj|2004h . 
Secondly, for numerical reasons the innermost bin is not re- 
solved in $, but averaged over all # € [0°, 90°]. 

In the spatial region with kinematical data the Coma 
galaxies offer different degrees of minor-axis anisotropy, 
from strongly tangential (GMP5279) to moderately radial 
(GMP3792). Towards the centre fig -> 0, while fi <p becomes 
negative (most likely due to the incomplete orbit sampling) . 
Going outward, many but not all galaxies exhibit a gradual 
change in dynamical structure, often in form of a minimum 
or maximum in fi. Around the last data point most models 
are isotropic. The most radial system is GMP3792. 

6.2 The equatorial plane 

Velocity anisotropy profiles in the equatorial plane are shown 
in Fig. [7] Note that unlike along the minor-axis axial sym- 
metry does not imply any relationship between fi v and fi$ 
at low latitudes. 
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Figure 6. Minor-axis anisotropy profiles. Top: meridional anisotropy bottom: azimuthal anisotropy solid: best-fit models 
(in the colour- version radial anisotropy is highlighted in red, tangential anisotropy in blue); shaded: 68 percent confidence region; 
dotted: region with kinematic data; arrows: r e g. 
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Figure 7. As Fig. [6] but for the major axis. 
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Figure 8. Short-axis anisotropy against minor-axis H4. The line 
in the upper panel shows a linear fit (quoted in the panel). 



Figure 9. As Fig. [8] but for the majc 



Meridional anisotropy. Contrasting the situation 
around the poles, almost no galaxy exhibits tangential 
anisotropy (3$ < 0. Apart from the peculiar object GMP3329 
(cf. Sec. 14- 1 p all galaxies have /3# > over the kinematically 
sampled radial range. The average (3-a turns out to be re- 
lated to the intrinsic flattening of the galaxies (Thomas et al. 
2007a, in preparation). Uncertainties on the intrinsic shape 
therefore propagate into uncertainties on /3#. As it has been 
stated in Sec. 13.11 in many cases it is not possible to distin- 
guish between different inclinations with high significance. 
Hence, intrinsic shapes are likewise poorly constrained and 
the uncertainties on /3# become large. A typical example is 
GMP3792: the best-fit inclination is i = 60° and requires 
a relatively flattened intrinsic configuration with large /3#. 
However, models at higher as well as lower i are within the 
68 percent confidence region. Consequently, the shaded area 
includes also models with different flattening and (3&. 

Azimuthal anisotropy. More diversity than in /3# 
is offered by azimuthal velocities. In GMP3792, for exam- 
ple, Ptp < suggests that the system may be composed 
of two flattened subsystems with low net angular momen- 
tum, causing large (^-motions. GMP3510, GMP3958 and 
GMP0144 are relatively isotropic (ay « a v ) over the kine- 
matically sampled spatial region. GMP5279, instead, offers 
A? ~ j3ip > 0, implying a r > a v and ay > a# over the region 
with data. 



6.3 Relations between anisotropy and observed 
kinematics 

The intrinsic short-axis velocity anisotropies are closely re- 
lated to the observed (local) H4. This can be taken from 
Fig. [5] where for each (projected) radius R with a measure- 
ment of H4, the local Hi{R) is plotted against the internal 
anisotropy P(r = R) at the same radius. Internal radii r 



have not been corrected for inclination since most models 
are edge-on (cf. Sec. I4.2|l . From the figure a tight correlation 
of (3$ with H4, follows (quoted in the plot): the smaller H4, 
the more tangentially anisotropic the model. A similar trend 
occurs between f3 v and H4 (lower panel). This reflects that 
A> ~ Pip around the symmetry axis (see above). 

Comparable trends between j3 and H4 hav e also 
been found in spherical mo dels (e.g. I Gerhard! Il993l : 
iMagorrian fc Ballantvnel l200ll ). The similarity between 
spherical models on the one hand and the polar region of 
axisymmetric models on the other might be connected to 
the fact that in both cases u v = o$. In other words, ef- 
fectively there is only one degree of freedom in the stellar 
anisotropy {f3 v = f3$) and, if the potential is fixed, there 
must be a close relationship between j3 and the shape of 
the LOSVD (as measured by H4). Experiments with spher- 
ical models indicate that the dependency of H4 o n the po- 
tential is weaker than its variation with /3 (e.g. I Gerhard] 
Il993l ; IMagorrian fc Ballantvnel 120011 ). If the same holds for 
axisymmetric potentials, then this would explain why along 
the polar axis of axisymmetric models ft depends in about 
the same way on H4 as in spherical models. Note, however, 
that our Schwarzschild models provide many more inter- 
nal degrees o f freedom than the smooth spherical models 
consid ered by I Gerhard! l|l99St ) and IMagorrian fc Ballantvnel 
(2001). This becomes apparent when the influence of regu- 
larisation on the fit is lowered and the scatter around the 
relation shown in Fig. [8] increases (cf. Sec. I9.2|l . 

In contrast to the polar region, no tight correlation be- 
tween H4 and velocity anisotropy is found around the equa- 
torial plane (cf. Fig. [9} . This holds especially for whereas 
there is a slight trend of fi v to increase with H4. For com- 
parison, a linear fit is shown in the lower panel. A detailed 
investigation of the orbital structure will be presented in 
another paper (Thomas et al. 2007a, in preparation). 
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Figure 10. Luminous matter phase-space densities. Top: prograde orbits are highlighted (blue); bottom: retrograde orbits 
highlighted (red); vertical dotted lines: boundaries of kinematical data; arrows: r c g. 
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7 PHASE-SPACE DISTRIBUTION FUNCTION 
OF THE STARS 

A more fundamental quantity related to a stellar dynami- 
cal system than its anisotropy is its phase-space distribu- 
tion function /. It describes the density of stars in phase- 
space and offers the most detailed and comprehensive view 
on its dynamical state. For stationary systems the DF 
is a function of the (isolating) integrals of motion and, 
thus, constant along indi vidual orbits (Jeans theorem; e.g. 
iBinnev fc TremaineHl987r ). To be considered in axisymmet- 
ric potentials are the energy E, angular momentum L z along 
the axis of symmetry and, in most astrophysically relevant 
potentials, the so-called third integral Iz- To be physically 
meaningful the DF has to obey the further condition that 
it is positive everywhere. In Schwarzschild models the con- 
stancy of the DF along orbits is explicitly taken into ac- 
count during the orbit integration. Its positive definiteness 
is guaranteed as long as the orbital weights Wi are positive 
(cf. equation [7]) . In other words, the very existence of our 
Schwarzschild models ensures that the luminous component 
of the model is stationary and physically meaningful (posi- 
tive density). 

A detailed investigation of the full dependency of the 
DF on all integrals of motion and its connection to stellar 
population properties will be the subject of another publi- 
cation (Thomas et al. 2007a, in preparation). Here we only 
consider some general properties of the DF. For this purpose 
it is convenient to define a mean orbital radius 

<r orb ), EE J2 ^ r >> ( 17 ) 

k * 

where T; is the total integration time of orbit i and r\ is its 
radius at time-step k (lasting Aif). In rough terms (r or b) can 
be interpreted as a measure of the orbital binding energy. 

Fig. 1101 surveys the DFs of all 17 Coma galaxies. Each 
dot represents the phase-space density of a single orbit. To 
roughly trace the angular-momentum dependency of the DF 
prograde orbits with L z > are highlighted in the top panel 
and retrograde orbits [L z < 0) are highlighted in the bottom 
one. 

The figure shows that with decreasing galaxy mass dif- 
ferences between prograde and retrograde orbits in phase- 
space become more significant. This partly reflects an in- 
creasing importance of rotation in lower mass galaxies of 
our sample. Often, the highest phase-space densities of pro- 
grade orbits are nearly constant over some radial region (e.g. 
GMP3958 around r e ff, GMP2440 inside r < r off ). In many, 
but not all, rotating galaxies the dominance of prograde 
orbits comes along with a strong depression of retrograde 
orbits. In the outer parts of GMP5975, for example, retro- 
grade orbits have phase-space densities up to 10 orders of 
magnitude smaller than prograde orbits. Such low-density 
orbits can actually b e regarded as being absent in the models 
l|Thomas et al.| [2005). Concerning the significance of this de- 
population it is interesting to note that in case of GMP5975 
it was originall y found in models ba sed on major and minor- 
axis data only dThomas et al.l i2005^ . but remains almost un- 
changed in our new models including additional kinematical 
information along a diagonal axis (cf. Sec. [2|. The depop- 
ulation of retrograde orbits cannot be a general modelling 
artifact since it does not appear in all rotating galaxies. A 



counter-example is the least-massive object, GMP3958: it 
rotates but does not show a strong depression of retrograde 
orbits in its outer parts. 

One galaxy of th e sample, GMP5568 , hosts a counter- 
rotating central disk l|Mehlert et alj|l998f ). which shows up 
by a dominance of retrograde orbits around (r or b) ~ 1 kpc 
and w/V « 10 M /pc 3 /(km/s) 3 in Fig. [TO] In most rotat- 
ing galaxies the majority of retrograde orbits follows ap- 
proximately a power-law like straight line (e.g. GMP3958, 
GMP2417). Retrograde orbits in GMP1990 follow such a 
power- law like distribution only outside r e s. Near the cen- 
tre of the galaxy retrograde orbits exhibit some excess den- 
sity, compared to a power-law extrapolation of the behaviour 
outside r e ff. This may reflect a faint inner counter-rotating 
sub-component, although unlike in GMP5568, prograde or- 
bits always dominate in GMP1990 (and the observed sense 
of rotation is the same at all radii). 

In some galaxy models orbital phase-space densities 
are more spread than in others: for example the DFs of 
GMP3329 and GMP5568 look particularly noisy. In case 
of GMP3329 its peculiar kinematics, already discussed in 
Sec. 15.11 may be responsible for the distorted phase-space 
distribution of orbits. 



8 PHASE-SPACE DISTRIBUTION FUNCTION 
OF DARK MATTER 

So far we have only considered the phase-space distribution 
function of the luminous component of our models. To en- 
sure that these models are physically meaningful we also 
need the dark halos to be supported by an everywhere posi- 
tive phase-space distribution function. Without the baryons 
present, the existence of DFs for our halo profiles is known. 
In case of NFW-halos it follows trivially from the fact that 
they arise in iV-body simulations a nd DFs for L OG-halos 
have been constructed explicitly by lEvansI (|l993T ). With a 
significant contribution of baryons (or any other component) 
to the overall gravitational potential, the existence of these 
DFs is no longer guaranteed, however. For example, if a 
cored halo (central logarithmic density slope 7 = 0) is em- 
bedded in a cuspy baryonic component (7 = — 1) and if the 
core radius exceeds a critical limit around rc > 3r e ff, then 
central phase-space densities become negative in isotropic or 
radia lly anisotropic systems l|Ciotti fc P ellegrini 1992; Ciottil 
Il999h. In co ntrast, a cuspy halo can always be supported 
l|Ciottilll996f ). Thus, the existence of a plausible halo DF for 
our LOG-halos, which often have core radii near or beyond 
the critical limit (cf. Tabs. [T] and [2]) is not obvious. The main 
goal of this section is to investigate whether we can find a 
positive definite DF for all our best-fit models, or whether 
the phase-space analysis rules out some of our halo profiles. 

8.1 Construction and existence of the halo 
distribution function 

Our modelling machinery allows to construct a DF for dark 
matter in an analogous way as for luminous matter: by solv- 
ing equation ([5]) with an orbit superposition. The only dif- 
ference to the calculation of the luminous matter orbit su- 
perposition is that now the dark matter density profile is 
used as the boundary condition and not the deprojected 
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Figure 11. As Fig, 1101 but for dark matter. No distinction between prograde and retrograde orbits is made (see text for details); 
large (blue) symbols: average phase-density of all orbits with (r OI ^,) < 0.lr e ff. 



light-profile v. In addition, since we lack of any kinematic 
information about the hypothetical halo constituents we set 
a = in equation (0 and maximise the entropy of the or- 
bits. For our goal of finding at least one positive definite DF 
this does not imply any loss of generality. 

Within the numerical resolution of our orbit models we 
find indeed orbit superpositions with positive orbital weights 
Wi > that allow to reconstruct the halo density profile in 
each case. The corresponding DFs are stationary by con- 
struction and positive everywhere. The fact that we even 
find positive definite DFs for those LOG-halos that are be- 
yond the above cited critical core-radius can have several 
re asons: our mo dels are slightly different from the ones used 
in lCiottil fl999) (different radial run in outer parts, baryonic 
component flattened in our case). In addition, our orbit su- 
perpositions can well produce tangential anisotropy, which 
helps to maintain a positive DF. Finally, our orbit models 
have a finite resolution. We cannot exclude that reconstruct- 
ing the halo density with higher resolution would force some 
orbital weights to become negative. 

8.2 Differences between NFW and LOG-halos 

Apart from the mere existence, there are significant differ- 
ences in the derived DFs, however. This can be taken from 
plots of the halo DFs in Fig. [Til NFW-halo DFs (GMP2921, 
GMP0144, GMP1176 and GMP5975) are monotonic with 
respect to (r or b) and regular. The high degree of regularity 
(compared to the corresponding luminous matter DFs) re- 



flects the maximisation of entropy, whereas noise in the stel- 
lar kinematics and sub-structuring of stars in phase-space 
tend to broaden the stellar DF (cf. Sec. [7J) . 

DFs of LOG-halos e xhibit a drop of central phase-space 
density, as predicted bv lCiottil l|l999t ). In GMP3329, where 
the halo is very concentrated, the drop is rather gentle. With 
increasing core-radius the drop becomes more substantial. 
In addition, the noise in the DF increases with increasing 
core radius. Such disturbances in the DF, even though we 
maximise the entropy, indicate that a fine-tuning of the or- 
bits is necessary for large core-radii to be supported by a 
positive-definite DF. This, and the non-monotonic depen- 
dency of orbital phase-space densities on (r or b) could imply 
that the corresponding DFs and, thus, also the spatial den- 
sity profiles are unstable. If this is indeed the case, then 
the phase-space analysis would provide a strong argument 
against large-cored halo profiles, independent from the kine- 
matic fits. Of course, the halo DFs shown in Fig. [TTJ are 
not unique, as stated above: the models have no access to 
the orbit distribution in the halo, apart from those con- 
straints coming from the shape of the density-profile alone. 
Details of the DFs in Fig. [TTJ are therefore physically mean- 
ingless. However, that the entropy maximisation does not 
yield smooth DFs for LOG-halos with large cores suggests 
that - independent from our ignorance about the details 
of the orbit distribution - smooth dark matter DFs in the 
corresponding baryonic potential wells are unlikely. 

In any case, a systematic stability analysis is out of the 
scope of this paper. What we can conclude here is, that 
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Figure 12. Best-fit Td yn versus regularisation parameter a. Dot- 
ted line: a = 0.02, the regularisation adopted for the best-fit 
models; shaded: one sigma confidence region for a = 0.02. 
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based on the kinematic fits and based on the mere existence 
of a positive definite halo DF, we cannot rule out one or the 
other halo profile with our orbit models. Shape and structure 
of LOG-halo DFs make them the less likely option, however. 



8.3 Central dark matter density 



According to Sec. 15.31 the central density of dark matter is 
often orders of magnitudes lower than the luminous mass 
density. This suggests that it is in many cases weakly con- 
strained. Could it be even lower than in LOG-halos? Accord- 
ing to the above phase-space analysis this seems unlikely, be- 
cause lower central densities would most likely augment the 
disturbances of the halo DF. Thus, the central spatial dark 
matter densities of our LOG-halos are likely lower limits to 
the true central dark matter densities. 



8.4 Central dark matter phase-space density 

From the dark matter DFs of Fig. [11] we have also calculated 
a mean central phase-space density 



(18) 



where the sums on the right hand side are intended to com- 
prise all orbits with (r or b) < 0.1r c ff. These central phase- 
space densities are flagged in Fig. [11] by the large symbols. 
Comparison with Fig. [10] shows that central dark matter 
phase-space densities are particularly low in systems with 
strong rotation. Exceptions are GMP1176 and GMP5975 
with their NFW halos. 

The uncertainty in the dark-halo DF related to our ig- 
norance about dark matter kinematics of course affects fh- 
According to our above discussion the drop in LOG-halo 
DFs seems a feature connected to the density profile, though, 
and we do not expect that reasonably isotropic or radially 
anisotropic LOG-halo DFs exist for which fh increases by 
orders of magnitude. Thus, although fh is subject to many 
uncertainties, it is likely good as an order of magnitude esti- 
mation for the central dark matter phase-space density con- 
nected with the mass decomposition made in Sec. 13.21 
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Figure 14. Best-fit circular velocity curves for different values 
of the regularisation parameter: a = 0.001 (blue, solid), a = 2.7 
(blue, dotted) and a = 0.02 (black, solid). 



9 REGULARISATION 

As it has been discussed in Sec. I3.4l the same regularisation 
a = 0.02 is adopted for all Coma galaxies. In the following 
we will discuss the dependency of our modelling results on 
the choice of a. 



9.1 The influence of regularisation on model 

masses 

Fig. [12] surveys the best-fit stellar mass-to-light ratios Td yn 
over the regularisation interval a G [10 ,3]. Two conclu- 
sions can be drawn from the figure. First, no systematic 
trend of Td yn with a is noticeable. In GMP1990, for exam- 
ple, Tdyn increases with a, while in GMP1750 it decreases. 
Second, in most of the sample galaxies the weight on reg- 
ularisation has barely any effect on Td yn (e.g. GMP5568, 
GMP0282, GMP0144, GMP5279, GMP2417, GMP3510, 
GMP1176, GMP5975, GMP2440). 

The best-fit dark matter fractions at three representa- 
tive radii are shown in Fig. [T3]as a function of a. As could 
have been expected, the dark matter fraction and Td yn are 
correlated: in most cases where Td yn , say, increases, the dark 
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Figure 15. Best-fit anisotropy at 0.1r e ff (black/solid), 0.5r e fj (red/short-dashed) and 1.0r e g- (blue/long-dashed); top: merid- 
ional anisotropy, bottom: azimuthal anisotropy; left: minor-axis, right: major-axis; vertical dotted lines: a = 0.02. 



matter fraction decreases (and vice versa). Since there is 
no systematic trend of Td yn with a it follows that there is 
also no systematic trend of the dark matter fraction with a. 
Moreover, the variation of dark matter fractions with a is 
within the quoted error budget of Fig. [S] 

Finally, Fig. 1141 shows circular velocity curves for three 
different values of a. The influence of a on the shape of the 
circular velocity curve is weak. Only in a few systems the 
general shape of the circular velocity curve changes with a 
(for example GMP3510). These changes occur mostly out- 
side the region covered by kinematic data, however. 

9.2 The influence of regularisation on model 
kinematics 

Now to the influence of a on the derived velocity 
anisotropies: the left panels of Fig. 1151 show best-fit merid- 
ional and azimuthal velocity anisotropies at three represen- 
tative radii as a function of a. The figures indicate that 
maximum entropy fits (a — > 0) yield isotropy along the 
minor-axis. Lowering the weight on regularisation gener- 
ally increases the anisotropy - the absolute value of f) — 
in the models, as could have been expected. There is no 
specific trend of /3 with a: some systems gain more tangen- 
tial anisotropy with increasing a (for example GMP5279), 
while others become more radial (for example GMP0756). 
In most cases the dependency of /3 on a is monotonic and 
/3(a) does not change sign. In other words, whether or not a 
galaxy model is tangentially or radially anisotropic does not 



depend on a. Only the exact degree of anisotropy changes 
with a. 

Major-axis velocity anisotropies are plotted on the right 
hand side of Fig. [TS] In contrast to the minor-axis case there 
is no trend of fig — » for a — *• 0. Variations of intrinsic 
velocity anisotropies with a are slightly weaker along the 
equator than they are along the minor-axis. Since the trend 
of j3 with a is again monotonic and sign preserving in most 
cases, the general property of a galaxy to be radially or 
tangentially anisotropic is insensitive to the particular choice 
of a. 

From the top-left panel of Fig. [TS] it is clear that the 
relation between (the a-dependent) anisotropy and (the a- 
independent) H4 described in Sec. HI] must change with dif- 
ferent amounts of regularisation in the models. This is illus- 
trated in Fig. 1161 which repeats the upper panel of Fig. [8] 
for two different values of a. The consequences of stronger 
regularisation are displayed in the top panel, while weaker 
regularisation leads to the distribution shown in the bottom 
panel. For comparison the linear fit from Fig. [5] is shown 
by the dashed line. Increasing the weight on regularisation 
makes the correlation tighter, but does not alter the slope. 
With less regularisation the scatter increases, because mod- 
els start to fit the noise in the data. The mean relation is in 
any case robust against different choices of a. Note, that a 
correlation between an intrinsic property (like f3$) and an 
observed one (like H4) cannot be the result of the entropy 
maximisation. In fact, for a = the meridional anisotropy 




Figure 16. As upper panel of Fig. \E\ but for different values of 
the regularisation parameter a (indicated in the panels). Dashed: 
linear fit for a = 0.02 (cf. Fig.©. 



Figure 17. Distribution of orbital weights for three differently 
regularised models (as indicated in the plot). All models are cal- 
culated in the best-fitting mass distribution of GMP5975. 



along the minor-axis vanishes ~ 0) and the relation 
breaks down. 

Concluding, the dynamical structure of the fits depends 
more strongly on the choice of a than the mass distribu- 
tion does (cf. Sec. I9.1|l . Thereby no clear trend of velocity 
anisotropies with a is noticeable. The monotonic behaviour 
of /3 with respect to a in most cases ensures that the general 
property of a model to be radially or tangentially anisotropic 
is independent of the choice of a. 

To give an example of how a influences the orbit dis- 
tribution Fig. [17] shows the histogram of orbital weights in 
the best-fit mass-model of GMP5975 for three different a. 
As can be seen the model at large a (weak regularisation) is 
dominated by a few orbits that carry almost the entire light. 
All other orbits are essentially depopulated in the model 
(only orbits with weights log to > —15 are included in the 
plot). The model at a = 0.02 is still relatively close to the 
maximum entropy distribution. 



10 DISCUSSION AND SUMMARY 

We have surveyed axisymmetric Schwarzschild models for 
a sample of 17 Coma early- type galaxies. The models are 
fitted to measurements of line-of-sight velocity distributions 
out to 1 — 4r e ff. Stellar mass-to-light ratios and dark halo 
parameters are determined for two parameterised halo fami- 
lies with different inner and outer density slopes. The models 
are regularised towards maximum entropy. 



10.1 Luminous and dark matter 

In each galaxy, models with dark matter fit better than mod- 
els without. A constant mass-to- light ratio is significantly 
ruled out in about half of the sample (eight galaxies where 



dark matter is detected on at least the 95 percent confidence 
level) . In four galaxies the case for dark matter is weak. The 
mass distribution in one of these systems (CMP 1990) is in 
fact consistent to follow the light. Five galaxies are inter- 
mediate cases where the formal evidence for dark matter 
is low, although fits with and without dark matter differ 
systematically in either their radial dispersion profiles or in 
their outermost LOSVDs. We believe that the low signal 
for dark matter in these systems is partly due to our very 
conservative treatment of the error bars. 

Our inferences about dark matter are based on the mass 
decomposition provided by equation JTJ. In particular, we 
have assumed that the stellar mass-to-light ratio is constant 
throughout the galaxy. In this context, GMP1990 is not sur- 
rounded by significant amounts of dark matter. However, 
before finally concluding upon dark matter in our galaxies a 
detailed comparison to independent estimates of the stellar 
mass is required (Thomas et al. 2007a, in preparation). For 
example, GMP1990 has a large mass-to-light ratio T = 10.0 
(J?c-band). If the actual stellar mass can only account for 
a fraction of it, then our result does not argue against dark 
matter in this galaxy, but merely implies that dark matter 
follows closely the light of the system. 

Constant mass-to-light ratios have been reported to 
be consistent with planetary nebulae kinematics in the 
outskirts of three round ish objects (spherical modelling; 
iRomanowskv et alj|2003l). Also some of the rou nd and non- 
rotating ellipticals of Kronawitter et all l|2000t ) are consis- 
tent with the mass distribution following the light distribu- 
tion. Many of these latter systems lack of kinematic data 
beyond r e g, however. In addition to the related uncertain- 
ties for the outer dark halo, spherical modelling of round 
galaxies generally suffers from the ambiguity related to the 
flattening along the line-of-sight. In this sense, GMP1990 is 
an interesting case, because its apparent flattening implies 
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a viewing- angle close to i = 90° and the kinematic data 
extend relatively far out (3r e ff). 

Best-fit dark matter halos are in 4 out of 17 cases of 
the NFW-type and in all other cases logarithmic. Differ- 
ences in the goodness-of-fit based on one or the other halo 
family are marginal in most cases. Central dark matter den- 
sities are at least one to two orders of magnitude lower 
than the corresponding mass densities in stars. Between 10 
and 50 percent of the mass inside the half-light radius r e g 
is formally dark in most Coma galaxies. These dark mat- 
ter fractions are in general agreement with earlier results 
of dark matter mode lling in round, non-rotating ellipticals 
jGerhard et alj|200ll) as well as with the analysis of col d gas 



kinemati cs dBertola et aL 19931; Qosterloo et alj 20021), hot 
halo gas llLoewenstein fc Whitelll999l ; lFukazawa et alj 20061; 
Humphrey et al.1 12006|) and st r ong lensmg stud i es l|Keetonl 
200ll ; lTreu fc Koopmanj2004h . lCappeUari et all I^OOrj 1 ) con- 
cluded for similar dark matter fractions in the SAURON- 
ellipticals (although in their models it is assumed that mass 
follows light). 

The combination of luminous and dark matter results 
in circular velocity curves of various shapes: some galaxies 
have outer decreasing w c i rc while others show an indication 
for a dip in v c i rc around ~ 10 kpc and a subsequent increase 
of Ucirc towards larger radii. We cannot easily quantify the 
significance of this dip. Its appearance close to the outermost 
kinematic radius is suspicious to reflect a modelling artifact. 
However, in contrast to the orbital structure, whose recon- 
struction becomes un certain around (and b e yond) the last 
kinem atic data point (|Krainovic et al ] |2005l : iThomas et al.l 
120051 ) , the mass reconstru ction in these regions is more ro- 
bust IjThomas et alj |2005). In addition, the dip does not ap- 
pear in all galaxies, suggesting that it is related to some 
observable property of the corresponding objects. In any 
case, it is interesting to note that similar dips are also indi- 
cated in temperature p rofiles of elliptical galaxy X-ray halos 



|Fukazawa et al.ll2006h and can also been seen in sp herical 
models of some round galaxies l|Gerhard et alj|200lf ). More 
extended kinematic data sets may in the future allow to bet- 
ter constrain the outer shape of the circular velocity curve. 

In rotating systems the circular velocity is fairly con- 
stant over the observationally sampled radial region (10 
percent fractional variation). Similarly flat circular veloc- 
ity curves have also been inferred from st ellar kinemat- 
ics of round systems (spherical mo delling; iGerhard et al.l 
I2001I : Im agorrian fc Ballantvne 200l [) an d from strong grav- 
itational lensing l|Koopmans et al.ll2006h . 

To the resolution of our orbit models all dark halos are 
supported by at least one positive definite phase-space dis- 
tribution function. In case of NFW-halos smooth DFs can 
be constructed, but for LOG-halos with large core radii even 
the maximisation of orbital entropy does not yield smooth 
DFs. It is not obvious whether the corresponding spatial 
density profiles are stable or not. Further modelling is re- 
quired to investigate whether phase-space arguments can be 
used to rule out logarithmic halos in elliptical galaxies. 



10.2 Kinematics 

With decreasing total mass the influence of rotation on the 
stellar phase-space distribution increases. In many galaxies 
rotation arises by an overpopulation of prograde orbits and a 



simultaneous underpopulation of retrograde orbits. At least 
one system lacks of the depopulation of retrograde orbits, 
proving that it is not a general artifact of our modelling 
approach. 

Some galaxies show strong tangential anisotropy along 
the minor-axis. This derives from low minor-axis H4- 
measurements, because observed H4 and modelled orbital 
anisotropy along the minor-axis turn out to be correlated. 
Slope and zero-point of this correlation are largely indepen- 
dent of regularisation, but - to some degree - stronger reg- 
ularisation tightens the relation. Such a relation between an 
intrinsic quantity on the one hand and an observed one on 
the other cannot originate from the entropy maximisation 
alone. 

Along the major-axis, a slight tendency of increasing 
(3 V with increasing H4 is noticeable, but with much larger 
scatter than along the minor-axis. Only one system shows 
indication of tangential anisotropy (f3$ < 0), all other galax- 
ies are isotropic or mildly radially anisotropic (/?# > 0, 
0if > 0). Radial anisotropy also a ppears characteristic for 
spherical models of round galaxies (|Kronawitter et ai1l2000l ; 
IGerhard et ai1l200ll ). A suppression of vertical energy (cor- 
responding to [5$ > 0) has recently bee n reported for 
SAURON ellipticals l|Cappellari et aljfeoOTt ). We plan a de- 
tailed investigation of the orbital structure of the Coma 
galaxies for the future. 

10.3 Regularisation 

Stellar mass-to-light ratios, dark matter fractions and the 
shape of circular velocity curves turn out to be robust 
against different choices of the regularisation parameter 
q. The strongest effect a has is on the reconstructed 
anisotropics: their absolute values tend to increase if the 
weight on regularisation constraints is lowered. At a fixed 
radius the dependency of anisotropy on a is mostly mono- 
tonic, such that the general quality of a galaxy to be radially 
or tangentially anisotropic, respectively, is independent of a. 
What changes instead is the actual amount of anisotropy. 

10.4 Outlook 

Detailed investigations of luminous and dark matter scaling 
relations, of stellar population properties and their connec- 
tion to the phase-space distribution of orbits are in prepa- 
ration. 
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APPENDIX A: DATA FITS 



Figs. I A 1 1 - 1 A 17l survey the fits to the photometric and kine- 
matical data for each galaxy (galaxies are arranged in order 
of decreasing total mass inside r e ff). 
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GMP3329 




Figure Al. Upper panel: Joint ground-based and HST photom- 
etry of GMP3329/NGC 4874. Lines: best-fit deprojection (red) 
and its edge-on reprojection (blue). Lower panel: stellar kinemat- 
ics along major (left /red) and minor axis (right/blue); filled and 
open circles refer to the two sides of the galaxy; dotted: best-fit 
model without dark matter. 

This paper has been typeset from a TfrjX/ JOTgX file prepared 
by the author. 
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Figure A2. As Fig. [All but for GMP2921/NGC 4889. 
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CMP55S3 GMP493& 




Figure A3. As Fig. |Al] but for GMP5568/NGC 4816. 
Green/third column: offset to major-axis r c ff/4; magenta/fourth 
column: offset to major-axis r c ff/20. 



Figure A4. As Fig. |XT] but for GMP4928/NGC 4839; 
green/third column: diagonal axis. 
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Figure A5. As Fig. El but for GMP0282/NGC 4952. 



Figure A6. As Fig.lATl but for GMP0144/NGC 4957. 
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Figure A7. As Fig. [Al] but for GMP5279/NGC 4827. 



CMF1990 




Figure A8. As Fig. ED but for GMP1990/IC 843; green/third 
column: offset to major-axis r e g/3. 
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Figure A9. As Fjg.lATl but for GMP1750/NGC 4926. 



Figure A10. As Fig. [All but for GMP3792/NGC 4860. 
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Figure All. As Fig. |Al] but for GMP0756/NGC 4944; 
green/third column: offset to major-axis r e g/2 (the two outer- 
most < —0.1 are omitted in the plot). 



Figure A12. As Fig. [Al] but for GMP2417/NGC 4908; 
green/third column: offset to major-axis r e g/2. 



CMF3510 
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Figure A13. As Fig. [All but for GMP3510/NGC 4869. 



Figure A14. As Fig. \M\ but for GMP1176/NGC 4931; 
green/third column: offset to major-axis r e g/3. 
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Figure A15. As Fig. ET] but for GMP5975/NGC 4807; Figure A16. As Fig.lAll but for GMP2440/IC 4045: green/third 

green/third column: diagonal axis. column: diagonal axis. 
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CMP3953 




Figure A17. As Fig. EH but for GMP3958/IC 3947. 



